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ABSTRACT 

A variety of observations now indicate that intergalactic helium was fully ionized by z ~ 3. The most 
recent measurements of the high-redshift quasar luminosity function imply that these sources had 
produced at least ~ 2.5 ionizing photons per helium atom by that time, consistent with a picture in 
which the known quasar population drives He II reionization. Here we describe the distribution of 
ionized and neutral helium gas during this era. Because the sources were rare and bright (with the 
photon budget dominated by quasars with luminosities L > L^,), random fluctuations in the quasar 
population determined the morphology of ionized gas when the global ionized fraction Xi was small, 
with the typical radius Rc of a He III bubble ~ 15-20 comoving Mpc. Only when Xi > 0.5 did the 
large-scale clustering of the quasars drive the characteristic size of ionized regions above this value. 
Still later, when Xi > 0.75, most ionizing photons were consumed by dense, recombining systems 
before they reached the edge of their source's ionized surroundings, halting the bubble growth when 
Rc ~ 35-40 Mpc. These phases are qualitatively similar to those in hydrogen reionization, but the 
rarity of the sources makes the early stochastic phase much more important. Moreover, the well- 
known characteristics of the z = 3 intergalactic medium allow a much more robust description of the 
late phase in which recombinations dominate. 

Subject headings: cosmology: theory - intergalactic medium - quasars 



1. INTRODUCTION 

For most of the Universe's history, the intergalac- 
tic medium (IGM) evolved rather slowly and smoothly. 
But there were two major exceptions: the reioniza- 
tion of hydrogen (at z > 6) and of helium (at 
z ~ 3). Recently, hydrogen reionization has re- 
ceived a great deal of attention in both the observa- 
tional and theoretical communities (see, e.g., reviews by 
Barkana fc Loebll2001l: ICiardi fc Ferrarall2005l : [Fan et all 
2006a; Fu rlanetto et al .l l2006bD . Helium reionization has 
received less attention, especially from theorists, but has 
actually been more accessible observationally and has 
important ramifications for quasar populations, galaxy 
formation, and the structure of the IGM itself. 

Because He II has an ionization potential of 54.4 eV, 
the relatively soft photons produced by (known popula- 
tions of) hot stars do not ionize it to any large degree 
(they can, on the other hand, singly ionize helium along 
with hydrogen). As a result, helium remained singly ion- 
ized until the quasar population built up in sufficient 
numbers: quasars have hard spectra, with luminosity 
densities L^, cx ly^^-^ (Tclfcr ct al. 2002) and can pro- 
duce enough photons to complete the reionizatio n pro- 
cess by z ~ 3 (ISokasian et al.ll2002t[Wvithe fc Loebll2003t 
iGleser et al.ll2005[ ). 

Observations do indeed point to helium reionization 
at z ~ 3. The strongest evidence comes from far- 
ultraviolet spectra of the He II Lya forest along the 
lines of sight to several bright quasars at z ~ 3: the 
observed wavelength of such an absorption system is 
A = 304(1 -I- z) A First suggested as a probe of the 
low-density IGM (jMiralda-Escudd 119931 : iGiroux et all 
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[T995I: iMiralda-Escude et allllQQfil : ICroft et al.lll997D . the 

He II forest has also proven to be a powerful probe 
of reioniza tion and of the background radiation field 
^Fardal et al. 1998). For our purposes, the most im- 
portant point is that the apparent He II optical depth 
decreases rapidly at z « 2.9, with a sprea d of Az » 
0.1 along differ ent li nes of sight ((ja kobscn e t al.l 119941: 



Davidsen et al."1996; Anderson et al. 1999; Heap et al. 
2000; Smct te et_al. 2002; Zheng ct al. 2004b; Shu U et all 
2004 iReimers et al.l 12004 120051: iFechner et al.l (2006^7 



This apparent transition is analogous to the rapidly- 
increasing H I optical depth ob served at z ^ 6 and usu - 
ally attributed to reionization (|Fan et al.ll200lL 120060 ). 
though see Becker et al. ( 200*3) for a different interpreta- 
tion. Although appropriate lines of sight are rare (they 
require a bright quasar with no intervening Lyman-limit 
absorbers that attenuate the far-UV fiux), recent large- 
scale surveys such as SDSS have opened up many more 
targets (e.g., Zheng et al. 2004a). 

We may thus be able to observe the helium reioniza- 
tion era in detail. Particularly interesting are the large 
opacity fluctuations ob s erved along severa l lines of sight 
jAnderson et al.' '199^; iHeap et al.l 120001 : ISmette et all 
[2002; Rcimers et al. 200^ At least one of these coin- 
cides with a nearby q uasar (Dob rzvcki fc Bechtoldl[l99l 
iJakobsen et al.l [2003( 1. This illustrates a crucial advan- 
tage of studying helium reionization over hydrogen: we 
understand the z ~ 3 Universe much better than the 
z ~ 6 Universe, in terms of the ionizing sources, the 
IGM, and other galaxies. We can thus provide much 
sharper tests of reionization models. 

There are several other indirect lines of evidence for 
helium reionization at z ^ 3, although in each case they 
are controversial. Helium reio nization shou l d at least 
double the IGM temperature. ISchave et al] (|2000( ) de- 
tected a sudden increase in the IGM temperature at 
z 3.3 by examining the Doppler widths of H I Lya 
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forest lines (see also ISchave et al.l Il999t iTheuns et alj 
l2002bl ): at about the same time, the equation of state 
of the IGM also appears to become nearly isot hermal, 
another indication of rec ent reionization (Schav e"et al.l 
120001 : [Ricotti et aD 12000 ') . However, temperature mea- 
surements via the Lya forest flux power spectrum show 
no evidence for any sudden change, although the er- 
rors are rather large (TZaldarriaga et al.ll2001l : TViel et al.1 
12004 iMcDonald et al ll2006[ ). Recent models of the he- 
lium reionization process predict sudd en jumps in the 
temperature as measured by the .Schave et al.l ()2000D 
method, although with somewhat smaller magnitudes, 
along with smoother evolution in the mean temperature, 
which may help to resolve th e controversy (|Gleser et al.l 
[2001 iFurlanet to fc Ohl[200l . 

Such a temperature jump should also decrease the re- 
combinat ion rate of hydrogen , thereby decreas i ng the H I 
opacity (|Theuns et al.ll2002aD . iBernardi et"al] (|2003f ) de- 
tected such a jum p at z ^ 3.2 in a l arge s ample of SDSS 
spectra; however, iMcDonald et alJ ()2006[ ) found no such 
feature with similar data an d precision. Most recently, 
iFaucher-Giguere et al.l (|2007D studied a sample of high- 
resolution Lya fore st spectra and found a feature nearly 
identical to that in IBernardi et al] (|2003f) . 

Finally, one would also expect the (average) metagalac- 
tic ionizing background to harden as helium is reion- 
ized, because the IGM would become transparent to 
high-energy photons. This should manifest itself in the 
He II /H I ratio (as pos sibly observed along one line of 
sight; iHeap et al]l2000( l and also in optical spectra. In 
particular, the ionization potentials of Si IV and C IV 
straddle that of He II, s o their rat i o should ev olve dur- 
ing hehum reionization. ISongailal ()1998L l2005| ) found a 
break in their ratio at z ~ 3; modeling of the ionizing 
background from optically thin and optically thick metal 
line systems also sh o ws a significant hardening a t z ~ 3 
(IVladilo et al.ll2003l : lAgafonova et al.ll20M l2007h . How- 
ever, other data of comparable quality show no evidenc e 
for rapid evolution ()Kim et al.ll2002l:lAguirre et al.ll2004D . 
This approach is made more difficult by the large fluctua- 
tions in the He Il/H I ratio eve n after helium reionization 
is complete (|Shull et al.ll2004l ). 

Despite this wealth of (often controversial) data, he- 
lium reionization has received relatively little theoreti- 
cal attention. There have been a few attempts to esti- 
mate the evolution of the globally-averaged He III frac- 
tion with redsh ift, given models for t he ionizing sources 
as input (e.g., IWvithe fc Loebl [2003h . There has also 
been one numerical simulation of helium reionization 
(jSokasian et al.l[2002l ) , which confirmed that quasars pro- 
vided enough photons to complete the process at z '--^ 3- 
4 and that the evolving opacity of the He II Lya forest 
was con sistent with the tai l end of reionization. Most 
recently, iGleser et al.l ()2005( ) used a semi-analytic Monte 
Carlo model to examine both the optical depth and tem- 
perature evolution in the context of helium reionization. 
However, there have been no efforts to understand the 
fundamental question of how the ionized and neutral gas 
are organized in the IGM (the "morphology" of reion- 
ization), which provides the overall paradigm in which 
observations must be interpreted. 

On the other hand, in the past several years hy- 
drogen reionization has received an enormous amount 
of theoretical attention, including both analytic mod- 



els and simulations (see iFurlanetto et al.ll2006bl for a re- 
cent summary). In particular, we now appreciate the 
importance of source clustering for reionization, which 
makes the process inhomogeneous on extremely large 
scales (> 10 Mpc) and substantial ly affects the inter- 
pretation of many obs ervables (iBarkana fc Loebl 120041 : 
IFurlanetto et al.ll2004bl hereafter FZH04). The purpose 
of this paper is to apply the simplest of these models to 
the helium reionization epoch as a first step toward un- 
derstanding the large-scale inhomogeneity of the process. 
In particular, we will describe three distinct stages in the 
evolution of the morphology, and we will emphasize the 
similarities and differe nces in the two eras. We base our 
models on FZH04 and IFurlanetto fc Oh! (|2005l hereafter 
FO05). 

This paper is organized as follows. We begin in ^j2] 
by considering simple reionization histories. We then 
describe stochastic and clustering-driven models for the 
bubble sizes in ^J3| and SjU respectively. We next consider 
the role of inhomogeneous recombinations in S|5l Finally, 
we conclude in ^ 

In our numerical calculations, we assume a cosmol- 
ogy with = 0.26, riA = 0.74, fit = 0.044, H = 
100ft. km s-i Mpc-i (with h = 0.74), n = 0.95, and 
as = 0-8, c onsistent with the most recent measurements 
(jSpergel et al. 2007). Unless otherwise specified, we use 
comoving units for all distances. 

2. THE REIONIZATION HISTORY 
2.1. An Empirical Model 

We first consider how the globally-averaged He III frac- 
tion, Xi, evolves with redshift. To begin, we will assume 
that quasars drive the reionization process. In that case, 
the key empirical input to our model is the quasar lumi- 
nosity function. We use the recent estimate of the bolo- 
metr ic luminosity function over a broad range of redshifts 
from lHopkins et all ((20073), who convert to a B-band lu- 
minosity function using the observed column density dis- 
tribution of quasars selected in the X-ray (thus account- 
ing for obscured sources). The distribution of _B-band 
luminosities L b (which we define to be i^L^, evaluated at 
4400 A, the center of the B-band'^) is shown in the left 
panel of Figure [T] for z = 3, 4, 5, and 6, from top to bot- 
tom. Here d$ /d log is is the number density of quasars 
per lo garithmic interval in lum inosity. 

The iHopkins et al.l (|2007cit ) luminosity functions are 
consistent with earlier estimates, but they do have sig- 
nificantly fiatter faint-end slopes (e.g. d^/dlogLs oc 
L° '^ for is < 10^2 Lq at z = 3). In the sample 
they used, the assumption of pure luminosity evolution 
no longer fit the data at high redshifts, w hich at the 
time showed a flattening at the faint end ([Hunt et al.l 
12004 ICristiani et"al1 12004D . More recent data show a 
somewhat steeper fa i nt end slope (oc L°-^ ) at z ^ 3 
(iFontanot et al.ll2007l : iBongiorno et al.ll2007l : isiana et all 
I2007D . more consistent with lower- redshift estimates. For 
simplicity we will use the IHopkins et al.l ()2007df ) function 
throughout, but we will also note where a steeper faint 
end slope affects our results. (It is not trivial to sim- 
ply change the faint end slope in the fit, because it is 
covariant with other parameters like L^,.) 

^ Note that we use a different luminosity convention than 
ISokasian etaLl ||200^ '). 
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Fig. 1. — Left panel: Quasar luminosity function from [Hopkins et al.l I l2007dl 'l. Right panel: Same, but weighted by Lg. In each panel, 
the solid, long-dashed, short-dashed, and dotted curves assume 2 = 3, 4, 5, and 6, respectively. Along the upper axis, we show the bubble 
size (in comoving units) corresponding to each luminosity, parameterized in terms of the fraction of ionizing photons absorbed inside the 
bubble (/abs), the quasar lifetime in units of 10^ yr (iy), and ^j^j^jj- 



The right-hand panel of Figure [T] shows d^/dLs 
weighted by luminosity, which is proportional to the 
number of photons produced per logarithmic unit of 
quasar luminosity. This panel explicitly pinpoints the 
quasars that contribute to the ionizing photon budget 
during reionization. Clearly, at z ^ 3, quasars with 
Lb ^ Li, Ri 10^^ Lq dominate. At higher redshifts, the 
luminosity function flattens at the bright end; hence, a 
much broader range of quasar sources contribute signifl- 
cantly. (However, note that the faint end converges even 
at z = 6, so quasars fainter than our assumed minimum 
quasar luminosity of 10^ contribute negligibly.) 

We must next convert from these B-band luminos- 
ity functions to the rate at which the quasars produce 
helium-ionizing photons with E > 54.4 eV: this requires 
a template for the spectral energy distribution of quasars. 
We use 



Li, oc 



-0.3 
-0.8 



2500 A < A < 4400 A 
1050 A < A < 2500 A 
A < 1050 A. 



(1) 



At A > 1050 A , this template agrees with 
that of iMadau et al.l (Il999( l: other templates (e.g., 
iSchirber k Bullock! I2003f disagree in detail but do not 
affect our conclusions, given the uncertainties. Most im- 
portant for us is t he far-ul traviolet spectral index a. At 
low redshifts, Telf er et al.l (j2002) find a wide variety of 
quasar spectral indices in the extreme ultraviolet, with 
a mean value of (a) « 1.6 and a standard deviation 
an ~ 0-8. This is slightly harder than the estimate of 
IZheng et all (|1998D . who found (a) « 1.8. We will use 
a = 1.6 for most of our calculations, but note that these 
uncertainties do affect our estimates of the absolute rate 
at which ionizing photons are produced. 
Given the template, the rate at which a quasar pro- 



duces ionizing photons is 

Lb ( T2&k 



iV,: = 0.0948- 



= 2.0 X lO'^'' s 



h V 1050 A/ 
Lb 



(2) 



(3) 



where in the second line we have assumed a — 1.6. (For 
a = 1.8, the prefactor becomes 1.4 x 10^^ s^^.) 

To develop intuition, it is convenient to assume a 
quasar lifetime iqso = lO^ty yr. In that case, a quasar 
can ionize a region of radius 



i?, » 14 



fabsh 



^HoII 



1012 L 



o 



1/3 



Mpc, 



(4) 



where /abs is the fraction of ionizing photons absorbed 
within the ionized bubble (as opposed to escaping into 
the surrounding IGM; see below) and the quasar lifetime 
is tqso — lO^h yr. We also assume that the bubble ex- 
pands into a uniformly ionized background with ionized 
fraction = l — x'^^n- Along the upper axes of FigurelU 
we show the bubble radius corresponding to each lumi- 
nosity. Clearly, we expect the typical bubbles around 
isolated quasars to be ~ 10-20 comoving Mpc across, 
although clustering will substantially increase the true 
sizes. 

We emphasize that Ri is considerably smaller than the 
nominal "proximity zone" around each quasar (i.e., the 
region in which its ionization rate exceeds that of any 
other quasar), even if the IGM is fully ionized. For ex- 
ample, the mean distance between L^, quasars at z = 3 

— 1/3 

is £ ~ Tiq ~ 110 Mpc ^ Ri. This is a reflection of the 
short lifetimes of quasars: [Ri/tf ~ t'jH{z = 3), so that 
over the age of the universe the entire quasar population 
can reionize the IGM. Most He III regions are not "ac- 
tive" but are instead fossils, which have been ionized by 
a now extinct quasar (see i)5.6p . 
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The ionized fraction Xi will evolve following 



dt 



ALf 



(5) 



where nfjo is the mean helium density, C = / {ne)"^ 
is the average clumping factor in the ionized IGM, and 
Au is the recombination rate per helium atom in gas at 
the mean density, 



An = aA{T)ne 



(6) 



We use the recom bination coefficients from 
IStorev fc HuWmI (fT995l) . ua = 2.18 x IQ-^^ g~i 
and as = 1.53 x IQ-^^ cm^ s'^ at T = 20,000 K. 
Unfortunately, the choice between case-A and case-B 
(which ignores recombi nations to the groun d state) 
is not an obvious one (jMiralda-Escudel ^OOl . Many 
studies use case-B, because any recombinations that 
regenerate an ionizing photon do not produce a net 
change in the ionized fraction. However, in a clumpy 
universe most recombinations actually occur near 
dense, mostly neutral systems (Lyman-limit systems for 
hydrogen- ionizing photons). Ionizing photons produced 
in these systems will be absorbed before they can escape 
into the low-density, ionized IGM for which we wish to 
compute the recombination rate. We therefore generally 
use the case-A rate below. Throughout this work, 
we will assume an IGM temperature T — 20,000 K, 
appropriate for photoionization with a relatively hard 
ionizing source. 

For gas at the mean density, the ratio of the helium 
recombination time to the Hubble time is 



i/-i(z) 



0.4 



l + z 



3/2 



aA{2 X 10* K) 



(7) 



where we have assumed matter domination. Because re- 
combinations are so fast, the sink term in equation ([5|) - 
and hence the clumping factor - is quite important. We 
will argue below that C will increase throughout reion- 
ization (FO05), but for the purposes of this simple cal- 
culation we will consider a range of constant clumping 
factors. 

FigurelSla shows ionization histories with C = 0, 1, and 
3 (dotted, solid, and dashed curves, respectively). In all 
cases, we use the iHopkin s e t al. (2007d) luminosity func- 
tions to compute the emissivity, as well as an ultraviolet 
spectral index a = 1.6. The ordinate, Q{z)^ is simply 
the time integral of equation ((51) ; it is the total number 
of ionizing photons produced per helium atom, minus 
the number of recombinations per atom, and can exceed 
unity once helium is fully reionized."' 

Given the observational evidence that helium was 
reionized at z ^ 3, it is reassuring to find that the 
iHopkins et al.l (|2007df ) luminosity function yields Q = 1 
at z ^ 3-3.6 in these models. There are a number of 
reasons why there may be a delay; for example, some of 
the ionizing photons have such high energies that they 
will not be absorbed in the diffuse IGM, our assumed 

* In a realistic model, C must increase rapidly as 5;^ — > 1, because 
then dense pockets of gas must begin to be ionized (FO05). This 
will of course keep < 1, but we do not attempt to enforce this 
limit. Thus, the very final stages of reionization may be more 
extended than depicted in our models. 




Fig. 2. — Helium reionization histories. (a): Using the 
IHopkins et al.l l|2007dl ) luminosity function, assuming C = 0, 1, 
and 3 (dotted, solid, and dashed curves, respectively), (b): His- 
tories normalized so that zjjc = 3. The thick curves (which all 
overlap) are identical to those in (a). The upper and lower dot- 
dashed curves follow /coii for rrimin = rrii and lOrrii, respectively; 



the thin dotted curves assume C, oc 



2/3 



with these same m,. 



spectrum may be too soft, and the clumping factor in- 
creases throughout reionization, especially when Xi ~ 1. 
Another uncertainty is in the faint-end slope, which 
more recent data shows to be somewhat steeper than as- 
sumed here (iFontanot et al.|[2007l : iBongiorno et al.ll2007l : 
iSiana et al.l 120071) . If we fix all the other parameters in 
the luminosity function, changing the faint end slope to 
have a constant shape oc L*^'^ increases the total lumi- 
nosity density by ~ 50%, 40%, and 35%, respectively, at 
z = 3, 4, and 5 (the difference occurs because the bright 
end slope flattens toward high redshifts) . Thus a steeper 
faint end will mean more photons, but it will not sub- 
stantially change the shape of the curves in Figure 

Figure shows reionization histories normalized so 
that Q{z — 3) = 1 (i.e., reionization completes at zhc — 
3). The thick curves, which nearly overlap, correspond 
to the histories in Figure Interestingly, the shape 
is nearly independent of the recombination rate, so long 
as C is constant and within this range. If we rewrite 
equation (O as dxi/dz, the sink term goes like ~ (1 -I- 
zy^^Xi{z); a fit to our results yields Xi{z) (1 -I- z)~^, 
much faster than the density and clumping evolution. 
Thus, to a good approximation the recombination rate 
simply tracks the rapidly-evolving ionized fraction and 
the shape of Q{z)/Q{z = 3) is nearly invariant. 

2.2. A Halo-based Interpretation? 

The empirical model described above gives our best 
handle on the history of Xi{z), but it does not provide 
enough information for the ionized bubble models that 
we examine next. These require a relation between the 
ionizing sources and the underlying large-scale density 

^ Another w ay to address t his qu estion would be to use the 
"Bright" fit of IHopkins et al.l Il2007dl ) , which fixes the faint end 
slope at L^'^. This also shows little effect on the shape of Q{z). 
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field. The clearest way to make this connection is to as- 
sign the quasars to dark matter halos. Although it is 
beyond the scope of this paper to construct a detailed 
model of quasar hosts, in this section we will use some 
general arguments to elucidate some aspects of their re- 
lationship. 

The thin dot-dashed curves in Figure [5^) use a different 
form for the emissivity, in which the quasar emissivity is 
assumed to trace /coii, the fraction of matter in collapsed 
halos. Thus the first term on the right hand side of equa- 
tion (O becomes 

^ d/coii 

Chc^, (8) 

where Chc is the number of helium-ionizing photons pro- 
duced per helium atom in these collapsed halos. For 
the minimum mass of a collapsed halo, we will take 
"Timin = CmTTii, whcre TTii corrcspouds to a virial tempera- 
ture Tvir = 2 X 10^ K, the approximate collapse threshold 
for halos in an ionized medium, and Cm is a redshift- 
independent constant accounting for the possibility that 
quasars only reside in massive halos. ^ For consistency 
with the excur sion set model outlined be low, we compute 
/coll using th e iPress fc SchechteH (1X9741 ) mass function. 
Although the lSheth fc TormenI mass function does 

match simulations better, the relation between halos and 
quasars is so uncertain that we regard the simpler form 
as reasonable (and more amenable to calculations). In 
gen eral, the abundance of high-mass objects is larger in 
the'S heth fc TormenI (fl999l) mass function, so for a given 
ionizing efficiency reionization would complete earlier. 

For quasar-like sources, the efficiency can be parame- 
terized as 



Cho — 7.3/ci 



/>54.4 150 eV £ /bh 
0.14 (i;ion) 0.05 10-5 



(9) 



where /esc is the fraction of helium-ionizing photons that 
escape the host galaxy of the quasar, />54.4 is the frac- 
tion of the quasar output emitted with E > 54.4 eV, 
(Eion) is the mean energy of these helium-ionizing pho- 
tons, e is the radiative efficiency of the quasar (so the 
total energy output is emeHC^, with ttibh the black hole 
mass), and /bh = ^n^Yi/^h is the fraction of the halo 
mass rrih inside the central black hole. To set these fidu- 
cial parameters, we have used the optical and ultrav iolet 
mean quasa r spec tra of IVanden Berk et al.l (|2001h and 
iTelfer et al.l ()2002[ ). respect ively, and the lo cal relation 
between ttt-bh and rrih from iFerrares^ ()2002l ) . 

The upper and lower dot-dashed curves in Figure dJ; 
show_models with mmin = "m-i and lOmi, respectively, 
and C — Q. Both are again normalized to complete reion- 
ization at z = 3, which requires Chc = 10 and 33 in the 
two cases. Interestingl y, these evolve significantly more 
slowly than the Hopki ns et al.l ()2007dl ) model, implying 
either that the black hole mass-halo mass relation must 
evolve with redshift or that quasars reside in more mas- 
sive halos than our simple halo model assumes (although, 
as we have emphasized, the faint end at z > 4 is quite 
uncertain, so the magnitude of the discrepancy is still 
unknown). A variety of other evidence points to a char- 
acteristic host mass of ~ 3 x 10^^-10^'^ Mq, suggesting 

^ Although in principle Cm could bo a function of redshift, we 
will assume for simplicity that it is a constant. 



that Cm ^ 10 is appropriate (e.g.. Figure 2 of iLidz et al.l 
[200fih . 

To develop some more intuition about the quasar-halo 
relation, it is useful to compare the quasar number den- 
sity to that of halos. At z = 3, the number density 
of quasars with L > 10^^'^ L© (roughly those domi- 
nating the photon budget, according to Figure [TJ)) is 
^ 4 X 10~^ Mpc~'^. Because each quasar is active 
for only a short time, the total number density of quasar 
hosts is nhost ~ nq[H{z)tQso]^^ 10"^^^^ Mpc"'^. De- 
pending on the assumed halo mass function, this requires 
^JT-min ^ 4-6 X 10^^ Mq ~ lOm^, roughly consistent with 
the reionization histories shown above. 

Of course, there are many lower luminosity quasars as 
well, although according to the empirical quasar luminos- 
ity function (Fig. [TJ)) they co ntribute a rela t ively s mall 
fraction of ionizing photons. iHopkins et al.l (|2005cl lallbl) 
have argued that lower-luminosity quasars are simply 
long, but relatively dim, phases in the same high-mass 
black holes that form high-luminosity quasars. Our fidu- 
cial model, with a high value of mmin, is meant to mimic 
the qualitative properties of this kind of quasar-host pre- 
scription. 

A second check is whether this simple picture repro- 
duces the o bserv ed clustering of high- redshift quasars. 
IShen et al.l (|2007l ) found that ^(r) « (r/ro)"'' with 7 = 
2.0±0.3and ro = 15.2±2.7/i~^ Mpc at z > 2.9, over the 
range 4 < {r/h~^ Mpc) < 150, for the relatively lumi- 
nous quasars observed with the Sloan Digital Sky Survey, 
with some evidence that the correlation length increases 
with redshift out to z = 4 even more strongly than one 
would expect from a volume-limited sample. In our cos- 
mology, this suggests a mean bias ~ 5.5 for the quasar 
halos, again corresponding to mmin ^ 6 x 10^"'^ M©. 

Despite this reassuring consistency, there is obviously 
a great deal more work to be done to model the relation- 
ship between quasars and their host halos. In particular, 
such a simple-minded picture does not quantitatively re- 
produce the quasar luminosity function as a function of 
redshift. In particular, simple one-to-one prescriptions 
tend to overpredict the abundance of faint quasars and 
produce a luminosity function that is too shallow at the 
bright end. However, it is reasonable to expect that a 
more sophisticated flavor of this class of models could: 
various modifications to the black hole-host halo relation 
allow one to simultaneously fit the lu minosity function 
and cluster i ng as a function of re dshift (jWvithe fc Loebl 
I2OO2I I2005t IHopkins et aHl2007bl |3). Next we illustrate 
some of the modifications req uired in such models. 

TheEHck ins et al.l (|2007af l model remedies this prob- 
lem by smoothly reducing the probability that both 
large and small halos host quasars, so that most quasar 
hosts are near a single characteristic mass; however, the 
width of this d istribution is not well-constrained (e.g., 
iLidz et an i2006). The motivation is that major mergers 
of large, gas-rich galaxies drive the formation of quasars, 
and these are mo st common in ho s ts jus t above a char- 
acteristic mass. iWvithe fc Loebl ([2001 use a similar 
merger picture to steepen the bright end of the lumi- 
nosity function. They also do not impose a minimum 
host mass and so produce a steeper faint end as well, 
more consistent with some of the recent data. 

Another problem is that, at z = 3, /coii(> lOm^) ~ 
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0.03, requiring a significantly larger ionizing efficiency 
than our fiducial choice in equation ([9]). Fortunately, 
there are several ways to alleviate this discrepancy. 
For example, many merger-driven models for quasar 
formation predict an increasing black hole mass frac- 
tion with redshift, ranging from m^Yil^h oc (1-1- 
z)°-^ (iHopkins et al.l |2 Q0ZbD to mBn/mh oc (1 -f zf^^ 
(jWvithe fc Loebl |2002f). (Both models appeal to quasar 
feedback to set the final relation but differ in how they 
relate the parameters of the central region surrounding 
the black hole to the global properties of the host halos). 
On the other hand, these scenarios would slightly worsen 
the discrepancy between the more slowly-evolving Xi{z) 
in these /cou-based models and the faster-evolving em- 
pirical calculation. 

Another possibility is that quasar luminosity 
scales more steeply than linearly with halo mass. 

iWvithe fc Lo"^ (jlOO^ argue that mBH « m'^J^ , over 
and above the redshift dependence. Because more 
massive halos form later, this will also help to accelerate 
the evolution in Xi. The thin dotted curves in Figure [2t 

show the difference: they assume ( oc m^f/^, with c„i — 1 
and 10 for the upper and lower curves. The acceleration 
is clear, and by weighting more massive galaxies more 
heavily we can indeed r each a shape v e ry clos e to the 
empirical one. However, 'Hopki ns et al.l (pOOTp) predict 
a much shallower relation with between black hole and 
halo mass (although still slightly superlinear) . The 
differing evolution of black hole mass with redshift in 
these models offsets the different black hole-halo mass 
relationships, allowing both to reproduce observations 
at the bright end. Both models provide at least approx- 
imate fits to the observed clustering (jWvithe fc Loebl 
[20051: IHopkins et al.l[20f)7c? ). 

Given that these more successful models appeal to 
mergers, is a simple halo-based interpretation really 
sufficient for our purposes? Instead, one might need 
to include only those halos that are growing rapidly 
(jCohn fc Chanel l2007f ) . Such a change would induce 
three changes in our /coii-based model. First, because 
major mergers are rare events, it would increase the 
amount of "stochasticity" in the quasar hosts relative to 
the halos. We will see in SjHjhow to describe this impor- 
tant aspect without appealing to either halos or mergers 
as hosts, so this aspect will not affect our results. Sec- 
ond, Q{z) could differ, since the halo merger rate and 
collapse rate evolve differently with time. Third, the dif- 
fering clustering of mergers could substantially modify 
the bubble distribution when it is driven by clustering. 
We examine these two effects below. 

The second effect is illustrated in Figure [3l which 
shows the rate at which galaxies accrete mass in merg- 
ers according to th e extended Press-Schechter model 
(|Lacev &: ColdfTOQSl) . More precisely, the curves show 



dm TO 



(Pp{M,m,z) 
dm dz ' 



(10) 



where d'^p{M,m,z)/dmdz is the rate at which halos of 
mass TO merge with a given halo of mass M per unit red- 
shift. We assume M > m and integrate over all mergers 
with mi < m < to„. We always take to„ — M (in order 
to avoid double-counting mergers); the upper and lower 
sets of curves take mi/M = 0.1 and 0.25, respectively, 




1012 10'= 
M (Mg) 



10"^ 



Fig. 3. — Fraction of original mass accreted by halos in mergers 
per unit redshift. The solid and dashed curves are at 2 = 3 and 4, 
respectively. All curves assume = M\ the upper and lower sets 
take rrt; = Af/10 and Af/4, respectively. 



showing two different possible thresholds for triggering a 
quasar phase. 

Figure |3] shows that the merger rate for halos of any 
mass increases with decreasing redshift, independent of 
the threshold. This means that merger models will 
steepen Q{z) relative to our fcow models, easing the ten- 
sion with the empirical model. 

The third effect is also illustrated in Figure |3l the vari- 
ation in the accreted mass fraction (normalized to the 
initial mass M) is surprisingly independent of M, ex- 
cept at masses well above the exponential cutoff in the 
mass function. For example, it varies by only ~ 5% over 
the entire W°-2 x lO^^ M© range at z = 3. Thus the 
range of halos that would contribute ionizing photons in 
a merger model is almost the same as in our fcoW models. 
The slight increase with mass for relatively small halos 



is much smaller than that in our ^ue 



2/3 

1 ' 



model. 



so such a picture would be bracketed by our two sets 
of models. The major difference with our fiducial /cou 
prescription would be to down-weig ht th e m ost ma ssive 
halos; this is w h y bot h the IHopkins et al.l (2007af ) and 
IWvithe fc Loebl |2002f ) steepen at the bright end. This 
is not too important to our picture, because the photon 
budget is dominated by halos near TOmin- 

Note that we use the extended Press-Schechter merger 
rates here; these have recently been shown to provide 
a relatively poor m atch to merger rates in the Mil- 
lenium simulati ons (iFakhouri fc Ma! l2007l ). However, 
iFakhouri fc Mai ([20070 find the same qualitative behav- 
ior in the simulations that we do; in particular, they 
show that the mean merger rate per halo depends very 
weakly on the halo mass (oc M° °^). Extended Press- 
Schechter overestimates the rates for equal mass merg- 
ers of 10"'^^ Mq halos by up to ~ 80 % at these redshifts 
(jFakhouri fc Mal[2?Ky7[ : iZhang et al.l [20081. However, the 
error decreases to ~ 40% for 10:1 mergers, and continues 
to decrease slowly for smaller mass halos. Scaling the 
merger rates in this way does not affect our conclusions. 
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because the dependence on mass is still gentle. 

In summary, matching Xi{z) suggests that quasars ap- 
pear only in massive halos (with Cm ^ 10) or that their 
luminosity scales steeply with host mass, although a de- 
tailed fit to the luminosity function requires much more 
machinery than we have introduced here. Given our un- 
certainties in the high-z quasar luminosity function and 
(especially) in the physics underlying the black hole-host 
relation, we will not attempt to construct a more de- 
tailed or accura t e mo del h ere; we simply note th at the 
IWvithe fc Loebl (|2002[ ) and lHopkins et al.l (|2007b[) mod- 
els are specific, and much more well-developed, cases 
of our more general picture. Instead we will exam- 
ine how this relation will affect the morphology of he- 
lium reionization. For instance, we explore the case 
where Chb oc rnff^, wh i ch rep roduces the main features 
of the IWvithe fc Loebl (|2002[ ) model, and a case where 
Cm = 10, roughly r e produc ing the characteristic mass of 
the iHopkins et al.l (|2007af ) model (although we do not 
directly appeal to mergers, as the original models did). 
Our qualitative conclusions will apply to any history, al- 
though some quantitative aspects will depend on the de- 
tails of the relation. 



2.3. A Stellar Contribution? 

While it is generally agreed that quasars drive helium 
reionization, it is worth noting that stellar systems can 
produce photons above the He II ionization threshold. 
In fact, the composite spe ctrum of z ^ 3 Ly man-break 
galaxies constructed by Shaplev et al.l (|2003f ) has a rea- 
sonably strong He II A1640 recombination line, even 
though it has little (if any) AGN contamination. The 
origin of this line is unclear. It is quite broad (FWHM~ 
1500 km s~^) and so is m ost likely atmospheric instead 
of nebular. However, Sha plev et al.l (2003) showed that 
known stellar populations, such as Wolf-Rayet stars, can- 
not simultaneously reproduce the He II line strength 
and the metallicity impl ied by other spectral features. 
iJimenez fc Haim"anl ()2006 ) have argued that the line may 
originate from massive Population HI stars that are able 
to form because of inefficient metal mixing. There have 
also been suggestions th at the He II A1640 line emission 
could arise coUisionally (|Yang et al.| [2006). 

In some of these scenarios, the He II A1640 luminosity 
may be accompanied by He II ionizing photons. This 
is the Balmer-a transition, so it must result from either 
recombinations following ionizations or collisional exci- 
tation to the n = 3 state. In the latter case, helium line 
cooling radiation peaks at T ~ 10^ K, where collisional 
ionization is weak; the gas tends to be optically thick to 
4 Ryd photons unless T > 10^ K ( Miniati et al.i 120041 ) . 
so there may not be much ionizing flux. But if the ori- 
gin is radiative, or if the gas is hot, these photons would 
be accompanied by a substantial ionizing flux. We now 
show that, if a substantial fraction of these photons es- 
cape into the IGM, they may play a significant role in 
helium reionization. 

Rather than work within the context of a specific 
model, we will instead simply use the IShaplev et al.l 
()2003f l spectrum to associate the He II A1640 luminos- 
ity, and hence the production rate of ionizing photons, 
to the star formation rate. (We could instead associate 
the luminosity to the total stellar mass, but it seems 



more likely that these high-energy photons are produced 
by young stars.) First, we note that the flux density at 
this line is /i, « 0. 35 /xJy, and the line has a n equivalent 
width Wx^ 2 A ^ imenez fc Haimanll2Q06D . Assuming 
that both kinds of photons result from radiative recom- 
binations, the line luminosity i s related to th e ionizing 
photon production rate Qi via (|Schaereill2003f) 



L ^ Cl(l - fcsc)Ql: 



(11) 



where /esc is the escape fraction of the photons and 
ci = 5.67 (6.40) X 10^12 gj.g ^j. g^g at 30, 000 (10, 000) K; 
we will assume the lower temperature for our estimates. 
From these , together with the average star formation 
rate in the IShaplev et all (|2003f ) sample (36 Mq yr^^), 
we find that the number of helium-ionizing photons pro- 
duced per helium atom in stars, rj, is 



70 



l-/c: 



2 A 



(12) 



The number of photons that escape into the IGM, per 
helium atom in the universe, is then 



]*{^) — Chc/coU — ??/csc/*/coll 



fa 



l-/e: 



/coll 

0.03 



(13) 

where Chc ^ measure of the ionizing efficiency per he- 
lium atom in galaxies and is the overall star formation 
efficiency. 

To estimate the factor /^/coii, we note that ~ 
10% of all bar yons are in stars at the present day 
(jFukugita fc Peeble s 2004) , and according to recent mea- 
su rements ^ 1/3 of th ese formed by z = 3 (e.g., see Fig. I 
of iPanter et al.l [20071 . and references therein). Thus at 
z = 3, /*/coii ~ 0.03. Comparing to Figure [2l this stellar 
component will be competitive with quasars if /osc ^0.5. 

The escape fraction is difficult to estimate, but it is 
possible (at least in principle) that it is larger than that 
for hydrogen (which is known to be at most a few percent; 
IShaplev et al.l[2006l ). For instance, the large line widths 
of observed Hell recombination lines 1500 km s~^) im- 
ply that He ionizing photons may arise from a different 
(or subset) stellar population - likely associated with 
stellar winds and outflows - from the general OB stel- 
lar population that produces hydrogen ionizing photons. 
Such a population with fast winds could possibly have a 
large fraction of "clear" sightlines. 

Thus we find that, even with the substantial helium 
recombination lines in moderate-redshift galaxies, stel- 
lar systems will probably not dominate the helium ion- 
ization budget, unless a large fraction (~ 0.5) of those 
photons can escape into the IGM. On the other hand, 
they may not be a completely negligible component, es- 
pecially since the resulting photons will probably be pro- 
duced quite close to the helium ionization edge (and thus 
interact in the nearby IGM) and because the stellar com- 
ponent evolves less rapidly than the quasar component 
(so may be more important at high redshifts). While it 
is certainly speculative to attribute helium reionization 
to this galactic emission, it could significantly impact the 
large observed hardness fluctuations in the ra diation field 
on sm all (~ 1 Mpc) scales after reionization (|Shull et al.l 
I2OOI . 

3. STOCHASTIC REIONIZATION 
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We now turn to modeling the morphology of helium 
reionization. We will begin by examining the limit in 
which the ionizing sources are relatively rare and bright, 
so that random fluctuations in their distribution deter- 
mine the ionization morphology. Of course, we do expect 
that quasars trace the underlying density field, so their 
locations are not truly random. But so long as the num- 
ber of sources per discrete ionized bubble^ is less than a 
few, random fluctuations in that number will dominate 
over clustering in setting the typical bubble size.^ This is 
different from hydrogen reionization, where the sources 
are smal l and numerous, so fluct uations can typically be 
ignored (jFurlanetto et al.l[2006al ). 

For a simple toy model, we begin by assuming that 
all the ionizing sources have an identical luminosity and 
that internal absorption within each ionized bubble is 
negligible. The first assumption is clearly not correct 
in detail, but at 2; = 3, ~ 50% of the ionizing photons 
are produced by quasars with L ~ 0.3-3-L*, so it is not 
a terrible approximation.^ We will examine the second 
assumption more closely below (see SJS|). These two as- 
sumptions then imply that the ionized bubbles are built 
of units with fixed volume Vi . 

We will further assume that the ionizing sources 
are P oisson-distributed with number density ngrc- 
iSheth fc Lcmson (1999) and Casas-Mira nda et all (|2002D 
found that the variance of halo counts in simulations of 
the z = universe is nearly Poissonian in regions ranging 
from voids to moderate overdensities, with the discrep- 
ancy relative to Poisson no more than a factor of two. 
Let us select regions of the universe with a smoothing 
window Vi ; we can compute the total ionized fraction by 
counting the smoothing windows that actually contain 
sources and multiplying by the number of sources within 
each such clump: 



x,^P{l\Vi) + 2P{2\V,) + 3P{3\Vi) + 

=S4 



fc=i 



(14) 



Here P{k\Vi) is the Poisson probability of finding k 
sources in each window; the prefactor k in the second 
line accounts for the extra volume ionized by each of 
these clumps,^*' and we have used ngrc 14 — Xi. (It is easy 
to verify that we recover the correct ionized fraction us- 
ing the power series expansion of e~^\) 

Similarly, we can compute the fraction of space in- 
side ionized bubbles containing at least two quasars by 

Sec ^4.1l below for a discussion of our discrete bubble approxi- 
mation. 

* Note that here we refer to the number of sources in a bubble, 
integrated over all time, not just the active quasars at any given 
instant. Most bubbles grow from multiple sources but will still 
have only one active quasar at any given time, because of their 
short duty cycle. We thus implicitly assume that the bubbles do 
not sign ificantly recombine between the episodic ionizations (but 
see ^5.61 1 . We also assume that a single halo does not host repeated 
quasar generations; in that case the bubbles would grow monoton- 
ically around these unusual halos and reach much larger sizes than 
estimated here. 

^ The fraction falls to ~ 35% with the steeper faint end 
slope of rec ent data (Font anot et al.| [2007: Bongiorno et al. 20071; 
ISiana et al.ll2007i ). Obviously the distribution of luminosities will 
be more important in this case. 

In other words, two nearby sources still ionize a volume 2Vi 
because of photon conservation. 




Fig. 4. — Fraction of the ionized volume filled by bubbles with 
at least N^i sources (top panel) and with exactly N^i sources 
(bottom panel), assuming that the they are Poisson-distributed. 
The dotted, long-dashed, solid, and dot-dashed curves assume 
Xi = 0.1, 0.3, 0.5, and 0.7, respectively. 



smoothing over windows 2Vi and including only those re- 
gions that are completely ionized. The same procedure 
works for bubbles with V > NVi, which must surround 
networks of N sources. Then 

H> N) = ^P(iVlV-) + ^;^PiN + m) + ... 



-E 



k=N ^ ' 

__ T{N)-{N-\)T{N-l,Nx,) 

' ^^^^ 

where the 1/A^ prefactor in each term appears because 
we are counting volumes in units of NVi. The fraction of 
space filled by bubbles with precisely N sources is then 
simply 

x,{N) = x,(>N)-x,{>N + l). (16) 

We show this distribution for several choices of Xi in 
Figured! The top and bottom panels show the cumula- 
tive and differential versions, respectively. Note that we 
normalize to the total Xi in each case, so the curves give 
the fractional contribution of each Nc\ to the total ion- 
ized volume. Obviously, the fraction inside small bubbles 
decreases throughout reionization, while the number of 
clumps grows fairly rapidly. We find that « 0.17%, 6.5%, 
and 24% of the ionized volume is filled by bubbles with 
at least five sources when Xi = 0.1, 0.3, and 0.5, respec- 
tively. 

Interestingly, for these randomly distributed sources 
there is no characteristic scale to the bubble size distri- 
bution; instead we find that the distributions approach 

One can imagine that, when Xi is sufficiently large, this 
smoothing procedure will not capture some configurations. How- 
ever, we will see below that the stochastic stage does not apply 
when Xi > 0.5 anyway. 
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power laws in the large N limit, with 



2x,{l-x,/2), 



(17) 



to a maximum absolute deviation of ^ 0.08 when Xi ^ 
0.35. We will see that this "scale-free" behavior is 
(in principle) different from reionization by clustered 
sources, which imprints a well-defined characteristic scale 
on the bubble distribution. However, in practice the bub- 
bles with one (or at most a few) ionizing sources domi- 
nate the distribution, so there will appear to be a charac- 
teristic scale of this magnitude in any realistic observed 
sample. This apparent scale is, however, illusory in an- 
other sense: we have ignored the luminosity distribu- 
tion of sources. In actuality, the size distribution will be 
broader as it must be averaged over the quasar luminos- 
ity function (though it will remain reasonably narrow be- 
cause Ri oc L^^^). This is a straightforward exercise for 
isolated quasars (see Fig.[T]), but it is difficult to include 
both a luminosity distribution and Poisson fluctuations, 
and we defer it to future, simulation-based work. 

Note that this simple model for stochastic bubbles does 
not make any assumptions about the quasar sources; it is 
based purely on the observed luminosity function. Thus 
it applies equally well to the halo-based model s we appeal 
to next and to merger-dri ven models (|Wvithe fc Loebl 



to next ana to mer ger-drr 
I2OOI iHopkins et al.lf2f)f)7ah 



4. IONIZED BUBBLES FROM CLUSTERED SOURCES 

We will now examine the bubble morphology generated 
by clustered ionizing sources. A simple way to see that 
clustering must be important is via the measured quasar 
correlation function. Since rp = (15.2 ± 2.7)h~^ Mpc 
for luminous quasars (|Shen et al.l 120071 ). the correlation 
length is comparable to the ionized bubble radius for an 
L, quasar, implying that clustering should certainly help 
large bubbles to build more quickly. Moreover, the clus- 
tering of halos near quasars certainly affects the local ion- 
izing background during hydrogen reioni zation ijYu fc Lul 
I2005t lAlvarez fc AbeH 120071 : iLidz et al.l [20071 ) . and it is 
reasonable to expect the same during helium reioniza- 
tion. 

For a more quantitative approach, we begin with the 
simple "photon-counting" model of FZH04. The basic 
idea is to compare the number of helium-ionizing pho- 
tons produced by the sources within a patch of the IGM 
to the number of helium atoms in that patch; wherever 
the former is greater, we have an ionized bubble. To 
begin, we assume a spatially homogeneous recombina- 
tion rate, with A^rcc the mean number of recombinations 
per ionized helium atom in the universe (but see 
If we further assume that every gravitationally-bound 
halo has a constant ionizing efficiency Cue, defined as in 
equation the criterion for a region with total mass 
M and fractional overdensity S to be ionized by sources 
contained inside it can then be written^^ 

CHc/cou(m„,i„, S\M) > (1 + iV,ec), (18) 

where /con(™min, S\M) is the collapse fraction in this re- 
gion. 

Note that, by using the cumulative number of recombinations 
per ionized atom, we are imphcitly assuming that the univ erse 
consists only of fully ionized and fully neutral regions; see il4.1l 
below. 



This condition p8|) can easily be modified to al- 
low the ionizing efficiency to vary across different ha- 
los (jFurlanetto et al.ll2006a[ ). If, for example, we allow 
Cne = CHe(mft), it can be written Cne/coU > 1 + Nrec, 
where 



Chc 



1 

/coll 



dmhCHc{mh)mhnh{mh), (19) 



and nii{mii) is the dark matter halo mass fun ction. 

We will use the Press fc Schechteil ()1974l ) model for 
nh{inh), although other mass functions give nearly iden- 
tical results when normalized to a constant average ion- 
ized fraction Xi, because our model depends only on the 
variation of the mass function with the large-scale over- 
density 5, w hich is more or less th e same for other mass 
functions (jFurlanetto et al.ll2006a[ ). Also, as discussed in 
' i2.2[ more detailed mo dels appeal to mergers r ather than 
halos to host quasars. ICohn fc Chand ()2007D presented 
a model, in the same spirit as FZH04, that follows hy- 
drogen reionization if sources are driven by mergers. In 
our case. Figure [3] shows that the merger rate is roughly 
independent of halo mass (except for extremely massive 
halos). We will see below that it is only the (luminosity- 
weighted) clustering that matters for our model, so this 
implies that basing the source population on halos and 
mergers will be nearly equivalent. Thus, we will defer 
a more detailed, merger-based treatment to the future 
and focus our initial model on the halo population, as in 
FZH04. 

FZII04 used condition (flSl) . toget her with t he excur- 
sion set formalism (Bond et a l. 1991; Lacev fc Col e 1995|) 
to write the mass function of ionized bubbles during 
reionization. In essence, for a_ region of mass M the con- 
dition Chc/coiK^I-^^) = (1 + -^rcc) (considered as a func- 
tion of the overdensity 5) replaces the usual spherical 
collapse criterion 5c = 1.69 for halo formation. This as- 
signs each excursion-set trajectory to the largest ionized 
bubble of which it is a part (thereby implicitly including 
all of the relevant neighbors). In order to write the solu- 
tion analytically, we replace the exact criterion for 5 with 
a linear fit in cr^ (m) , the smoothed linear-theory variance 
in the density field. The resulting mass distribution of 
ionized bubbles is 



?T.b(m, z) 



diner 



dlni 



Bo{z) 
a{m) 



exp 



20-2 (m) 



(20) 

where the excursion set barrier is w B{m^ z) — Bo{z) + 
Bi{z)a'^{'m) and p is the mean comoving mass density. 
Note that the mean recombination rate has no effect on 
the result (at a fixed Xi), because it is completely de- 
generate with the ionizing efficiency. (Inhomogeneous 
recombinations, on the other hand, have a substantial 
effect; see Sj5] below.) 

Thus the original FZII04 model is easily portable to the 
case of helium reionization, requiring only a re-definition 
of the ionizing efficiency. Figure [5] shows some example 
distributions. For our fiducial model (solid curves), we 
assume that rTimin = lOmi and that ^hc is independent 
of halo mass, which gives an ionization hist o ry clos e to 
(but slightly slower than) the IHopkins et all ([20073) lu- 
minosity function (Fig. [2|). Panels (a) and (b) contrast 

this with models with Cho c>c m^' and with TOmin = rrii, 
respectively. 
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Fig. 5. — Bubble size distributions at 2 = 4. The solid curves 

'Ttmin = lOrrii. In (a), the dashed curves take rrimin = lOrrii with (^hc 
ionizing efficiency. In each set, the curves have Xi = 0.2, 0.4, 0.6, and 




R (Mpc) 

in both panels assume (^He is independent of halo mass and take 

oc m^yl^ ■ In (b), the dashed curves take m^\j^ = vrn with a constant 
0.8, from left to right. 



For our fiducial model, the distributions peak at 
i? - 4, 12, 30, and 60 Mpc at Xi = 0.2, 0.4, 0.6, and 
0.8, respectively. The bubbles are somewhat larger if 
^He is an increasing function of mass, and somewhat 
smaller if lower-mass halos contribute. This is be- 
cause the bubble sizes are primarily a function of the 
luminosity- weighted clu stering of the ionizing sources 
(jFurlanctto et al." 2006a). 

Figure [5] contrasts the distributions at z = A and z = 3. 
As with hydrogen reionization, the dependence on red- 
shift is extremely weak: the bubble sizes are sensitive to 
the integrated bias of the sources, which does not evolve 
rapidly in this regime. They are slightly smaller at the 
lower redshift, primarily because mt is a decreasing func- 
tion of redshift. 

Of course, we have already seen that for reasonable 
quasar lifetimes the ionized region around a single 
quasar at 2; = 3 has i? ~ 14 Mpc (they are somewhat 
smaller at larger redshifts because decreases); our typ- 
ical bubbles do not reach such large sizes until Xi ~ 0.4. 
It may at first seem that the smaller bubbles are just 
created by dimmer quasars that produce relatively few 
ionizing photons. However, there is actually a fairly fun- 
damental inconsistency here. 

The problem occurs because the clustering model as- 
sumes that the ionizing bubbles precisely trace the un- 
derlying density field and does not self-consistently ac- 
count for the discrete sources: in other words, it does not 
properly account for their stochastic fluctuations, and it 
forces high-density regions to be ionized even if they may 
not have a high-mass halo. This tends to force many 
small, high-density regions to be ionized, rather than the 
much larger environments of the real halos. As a re- 

In other words, these regions' average collapse fractions are 
large enough to ionize themselves, but not enough to contain an 
actual halo. For example, consider a region that is 20% the size of 
a typical ionized bubble, and that has a 20% chance of containing 
a halo with m > m^-m. According to the density driven model. 
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Fig. 6. — Bubble size distributions at 2 = 4 (solid curves) and 
2 = 3 (dashed curves). Both sets assume our fiducial model. In 
each case, the curves have Xi = 0.2, 0.4, 0.6, and 0.8, from left to 
right. 



suit, we must turn to the stochastic model of ^3] early in 
reionization. 

This can be seen most clearly by examining the char- 
acteristic bubble size Rc (defined as the peak of Rmrib) 
to an analogous quantity in the stochastic reionization 
model. As we have seen the size distribution in the lat- 
ter case is nearly a power law: there is no well-defined 
peak. We therefore compare Rc to the effective number 
of sources Nc for which V{> Nc) — Xi/2 (i.e., half of the 

every one of these regions would be ionized, even though in reality 
only one in five would be. But that one region would also ionize 
its surroundings, making the total ionized volumes equal in each 
picture. 
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Fig. 7. — Characteristic bubble sizes as a function of 5;^ at 2: = 3. 
The solid curve assumes i^He is independent of halo mass and takes 
™min = lOmi. The dotted curve takes m^^^ = nii. The short- 
2/3 

dashed curve takes i^ue ^ ™^ • The long-dashod curve is for our 
stochastic model (see text for a precise definition). 



ionized volume is in bubbles smaller than - or larger than 
- this value). For these purposes, we treat iVc as a contin- 
uous, rather than discrete, variable. We then convert to 

1/3 

an effective radius by assuming that Rc = 147Vc ' Mpc, 
where the proportionality constant comes from the peak 
of the z — i luminosity density in Figure [Ij) (see eq. 14]) 
and assumes /abs = Ij ^7 = !> and x^^y^ — 1. 

Figure [7] shows how these scales evolve in several mod- 
els. The solid curve shows our fiducial clustering model, 
the short-dashed curve assumes Chg oc to^^^, and the dot- 
ted curve assumes mmin = "i^. The long-dashed curve 
shows the stochastic model, with defined as above. 
As with hydrogen reionization, the characteristic size of 
the clustering-driven bubbles increases rapidly through- 
out reionization, surpassing 15 Mpc by Xi ^ 0.5 and 
eventually reaching Rc > 100 Mpc. At a given ionized 
fraction, R^ increases with the underlying bias of the 
sources. 

In comparison, R^ defined from the stochastic model 
increases much more slowly: by only a factor ~ 2 from 
Xi = 0.01 to Xi — 0.75. Thus, only a few sources con- 
tribute to each "typical" bubble, regardless of the ion- 
ized fraction. The density-driven model predicts much 
smaller bubbles early on, with only a small fraction of 
the Universe contained in bubbles above the stochastic 
scale. However, as Xi increases Rc increases much more 
rapidly for the clustering model than for the stochastic 
model, eventually overtaking it. During the latter stages 
of reionization, the stochastic component will become 
negligible compared to the typical ^ 100 Mpc bubbles 
that emerge; by this point, fluctuations in the clustering 
networks are no longer important. 

Of course, we have built our stochastic model from 
bubbles with a fixed size, appropriate for Lj, qu asars. 
This is reasonable with the iHopkins" et all (|2007dr ) lumi- 
nosity function, which is rather strongly centered around 
Li, at high redshifts. With the steeper faint end slope in- 



dicated by more recent data, fainter quasars contribute 
somewhat more of the luminosity density. Thus the 
stochastic distribution will be built of somewhat smaller 
units, and we would expect clustering to become impor- 
tant somewhat earlier. 

This suggests that we can approximately account for 
the stochastic component by assigning all of the volume 
nominally in small bubbles to the characteristic scale of 
the stochastic component. Of course, that distribution 
is relatively uninteresting (depending only on the mean 
number density and luminosity of sources), so it is only 
once the typical bubble grows beyond the stochastic scale 
that we can study the astrophysics of the sources via their 
effects on the bubbles. For example, we have seen that 
the bias of the quasar host population affects the bubble 
sizes. Alternatively, if there is a significant population 
of small bubbles early in the reionization process, this 
raises the interesting possibility that there is an unre- 
solved population (such as miniquasars or galaxies) that 
contributes to helium reionization. 

4.1. Bubbles of Ionized Helium? 

In the previous two sections, we have treated the ion- 
ized zones around quasars as having sharp, well-defined 
edges. This is necessary for the excursion set model of 
FZII04, which implicitly assumes that gas is either fully 
ionized or fully neutral. It is an excellent assumption for 
hydrogen reionization; the mean free path of a hydrogen- 
ionizing photon is Am w 0.002([l + z]/10)^(zy/fc'Hi)^ Mpc, 
where i^hi corresponds to the ionization edge of neutral 
hydrogen. The stellar sources generally assumed to be re- 
sponsible for hydrogen reionization have rather soft spec- 
tra, so (v) « z^Hi- Thus, once an ionizing photon hits the 
edge of an ionized bubble, it is quickly absorbed. The 
transition regions between ionized and neutral gas have 
negligible thickness, so the two-phase approximation is 
an excellent one. 

In the case of helium, several factors make this assump- 
tion problematic. First, helium is less common and has 
a smaller photo-ionization cross-section. Second, helium 
reionization occurs at 2; ^ 3, when the universe is much 
more dilute. As a result, in a homogeneous medium 

AHo = 0.66f-^) (^) Mpc, (21) 

where /it'Hcii = 54.4 eV. Third, and most importantly, 
the quasars thought to be responsible for helium reion- 
ization have rather hard spectra. At the mean spectral 
index, a = 1.6, only half the photons have v < 1.5!/hcII 
and hence have mean free paths Ahb < ^1/2 ~ 2.3[4/(l -|- 
z)]^ Mpc. Thus we expect the typical "transition region" 
around each ionized bubble to be a few comoving Mpc 
thick. Fortunately, we have seen that the bubbles exceed 
this size throughout most of reionization (see eq. H]), so 
our two-phase approximation should be reasonable. We 
also note that the o nly published full nu merical simula- 
tion of reionization (jSokasian et al.ll2002t l also made the 
two-phase approximation by assuming that all photons 
lay at the ionization edge. 

The hard photons that do escape to large distances 
well outside the transition region form a more uniform, 
low-level ionizing background. Bubble expansion in a 
uniformly ionized medium can easily be incorporated into 
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the FZH04 model by chang ing the ionization criteri on to 
Cho/co11 > ShcIi(1 + ^roc) ijFiirlanetto et al.l l2004a^ : the 
net result is to make the bubbles somewhat larger for a 
given total ionized fraction than they would otherwise 
be, but the qualitative behavior is unchanged. 

There is one additional complication: the IGM is, of 
course, clumpy and not uniform. We will next describe 
how this affects the bubble distribution. 

5. INHOMOGENEOUS RECOMBINATIONS AND THE 
BUBBLE DISTRIBUTION 

To this point, we have ignored the effects of recombi- 
nations on the bubble distribution (treating them as a 
spatially uniform photon sink). FO05 showed that inho- 
mogeneous recombinations become increasingly impor- 
tant during hydrogen reionization as denser and denser 
regions are ionized. They constructed a second excur- 
sion set barrier that corresponded to the (instantaneous) 
ionization rate balancing the recombination rate inside 
the bubbles (with the latter set by the requirement that 
ionizing photons be able to reach the edge of the bubble, 
which sets the minimum density of neutral regions via 
mean free path arguments). This requirement becomes 
important in the late stages of reionization and effec- 
tively caps the bubble growth; qualitatively, this maxi- 
mum size corresponds to the mean free path between op- 
tically thick Lyman-limit systems in the early universe. 
The transition to the recombination-limited regime en- 
forces a smooth match from the highly inhomogeneous 
bubble-dominated morphology characteristic of reioniza- 
tion to the more uniform radiation background of the 
post-reionization universe. 

The FO05 excursion set barrier implicitly assumes 
that the rate at which ionizing sources produce pho- 
tons is steady. This is an excellent approximation for 
the small sources expected to dominate the hydrogen 
era, especially if their st ar formation timescales are long 
(|Furlanetto et al.ll2006a^ . However, helium reionization 
is very different: as we have seen, bubbles are typically 
produced by only a few sources, and the duty cycle is 
very short: tQsoH{z) 0m3t7[{l + z)/4:f^^. Thus even 
the largest bubbles are likely to contain only one (or at 
most a few) active quasars at any given time. In these 
bubbles, the emissivity fluctuates strongly, at zero most 
of the time but far above the mean value during active 
phases. Thus, we really only need to consider IGM re- 
combinations around isolated quasars. 

To that end, consider a gas element in the IGM with 
overdensity A ~ p/p. In the highly-ionized limit, ioniza- 
tion equilibrium demands that 



r 



:3.8 X lO^'^AT, 



-0.7 



r 



,(22) 



where we have used the case-A recombination rate, T = 
10^T4 K is the IGM temperature, and the ionization rate 
is 

<Ji, 
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r = (i + z)2 / 
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Fig. 8. — Density required to produce an absorption system 
optically thick to photons at the helium ionization edge, as a func- 
tion of distance from the central quasar. The dotted, dot-dashed, 
short-dashed, long-dashed, and solid curves assume log(LB/L0) = 
9, 10, 11, 12, and 13, respectively. Along each curve, the solid tri- 
angle marks the point at which R = A(Ai) according to the MHROO 
model. We take z = 3. 



Here cr^ = 1.91 x 10~i*(fHeii/i^)^ cm^ is the photoioniza- 
tion cross section for He II, R is the comoving distance 
from the quasar, and we have used the mean quasar spec- 
tr um (a = 1.6). 

iSchavd (|2001j ) has shown that H I Lya forest ab- 
sorbers can be accurately modeled by assuming their 
physical scale to be comparable to the Jeans length 
Lj ~ Cs/{Gpy^^ . We will use the same approximation 
for helium, so that an absorber with overdensity has 
column density A^hoII ~ LjXiUnc- Then 



TVhott « 1.8 X lO^^A^/^T, 



(25) 

A system will become optically thick when T,y — 
Ciz-^Hoii = 1- At the ionization edge, this requires an 
overdensity A^ 



A,; « 64r; 



2/15 



2/3 /Tvj \ 4/3 /, , \ -5/3 



Mpc 



(26) 

Note that this definition of A^ uses the photoionization 
cross section at the ionization edge; with the hard spec- 
tra typical of quasars, many of the photons have much 
higher energies and so can penetrate even denser systems 
than implied by this simple model (see M.ll above). To 
account for this, one could define a frequency-dependent 
= Aiiy I VYio)"^ ■ On the other hand, many high- 
energy photons will not even interact at the edge of the 
bubble but instead join the uniform background. 

Figure [8] shows this overdensity A^ as a function 
of distance for quasars of several different luminosities 
(log[_LB/LQ] = 9, 10, 11, 12, and 13, from top to bot- 
tom). The faintest quasars have A^ < 5 outside of a 
few Mpc, but typical sources are able to ionize even ex- 
tremely dense gas up to rather large distances. 
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To estimate the effects of these opticaUy thick systems 
on the ionized bubbles, we will use the same fundamen- 
tal approach as FO05 but take advantage of the much 
better-known IGM at 2; ~ 3. In particular, if a region 
with A > Ai appears between a quasar and the edge 
of its host bubble, then most of its photons will be con- 
sumed overcoming recombinations rather than ionizing 
new material. We therefore wish to compare our bubble 
sizes to the mean free path between IGM patches with 
A = A^. To obtain the latter, we need a model for the 
density structure of the IGM. 

5.1. The IGM Density Distribution 

iMiralda-Escude et all (|2000l hereafter MHROO) 
present a model for the volume-averaged density dis- 
tribution of the IGM gas, Py(A), based on a fit to 
simulations at z '--^ 2-4: 



Py(A)dA = yloA-'3exp 



(A 



-2/3 



2(25o/3)2 



dA. (27) 



Intuitively, the underlying Gaussian density fluctuations 
are modified through nonlinear void growth and a power 
law tail at large A. Here Sq essentially represents the 
variance of density fluctuations smoothed on the Jeans 
scale for an ionized medium; thus Sq oc {1 + z)^^ at 
high redshifts. The power-law exponent f3 determines 
the behavior at large densities; for isothermal spheres, 
it is /3 = 2.5. We use the fitted values from MHROO 
for these parameters; the remaining constants {Aq and 
Co) can be set by demanding proper mass and volume 
normalization. 

Note that MHROO base their distribution on a simu- 
lation that does not incorporate all the physics of he- 
lium reionization, and in particular it does not follow 
the inhomogeneous heating during that event. In reality, 
the clumping will decrease durin g helium reionization be - 
cause of photoheating (see, e.g.. lFurlanetto fc Ohll2007l ). 
which will modify the distribution of densities as well. 
This is probably not too important for our purposes, be- 
cause the IGM takes a fairly long time to adjust to the 
new Jeans smoothing scale (jGnedin fc Huil [l998h . and 
the increase is much more modest than during helium 
reionization. 

MHROO also offer a prescription for determining Ai, 
the mean free path of ionizing photons. In their model, 
it equals the mean distance between clumps with A > A^ 
along a random line of sight, which is approximately 



A, = Ao[l-Py(A,)] 



-2/3 



(28) 



Here Fv{Ai) is the fraction of volume with A < A^ and 
Ao is a (redshift-dependent) normalization factor. For- 
mally, this expression is valid only if the number density 
and shape (though not total cross section) of absorbers 
is independent of A^. This is obviously not true in de- 
tail for the cosmic web. However, MHROO found that it 
provided a good fit to numerical simulations at z = 2-4 
if we set \qH(z) = 60 km s^^ (in physical units). 

5.2. The Recombination Limit 

For simplicity, we will assume that gas with A < A^ is 
highly ionized while all gas with A > A^ remains neutral, 
so that Ai is an accurate approximation to the mean free 
path. This provides a reasonable description of shielding 



in dense regions, if those regions can be considered to be 
isolated clumps in which the density increases inwards 
(see the Appendix to FO05 for more detail). In that 
case, the radiation field will ionize the outskirts of the 
cloud until r « 1. Because of the density gradient (in- 
creasing inwards) , this skin corresponds to our threshold 
Ai. In detail, equation ignores two significant (but 
opposing) effects. First, it probably overestimates Ai by 
up to a factor ~ 2 because of accumulated photoelectric 
absorption by low column density systems (FO05). But 
second, it ignores higher energy photons, which are able 
to travel farther. We will assume that these two effects 
roughly cancel. 

To obtain i?max, the maximum bubble size in the pres- 
ence of recombinations, we simply find the point at which 
the mean free path Ai(Ai) equals the distance from the 
quasar. The filled triangles in Figure [5] show i?inax for 
several quasar luminosities. It ranges from ~ 5-100 Mpc 
over this luminosity range. For the faintest quasars, 
Ai < 5; in these cases, the assumption of a sharp transi- 
tion in the neutral fraction at Ai is probably not a good 
one, because the gas is so close to the mean density. How- 
ever, at L > 10" Lq, a, > 8 and Fv{Ai) < 0.03, so the 
two-phase approximation is probably reasonable. 

There is one important caveat about this method: the 
MHROO procedure computes the mean free path at an 
arbitrary point in the IGM, whereas we are interested in 
the mean free path as seen from the qua sar, which most 
likely sits in an overdense regions (see a lso lYu fc Lull2005l : 
lAlvarez fc Abeill2007HLidz et al.ll2007f ). Thus we proba- 
bly overestimate the mean free path. However, even at 
z ~ 6 (where these massive halos are much more rare) 
the environments typically approach the mean density 
within < 20 Mpc of the quasar, so this should not be a 
huge effect. 

Figure [9] shows i?max as a function of luminosity for 
quasars at 2: = 3 and 4 (solid and dashed curves, re- 
spectively). The filled triangles indicate the luminosity- 
weighted mean i?max across the entire quasar popula- 
tion. The dotted curve shows Ri from equation ([4]), the 
maximum size of the ionized bubble surrounding an iso- 
lated source (assuming t-j = \, Xngji — 1, and /abs = !)• 
We find (i?max) ~ 35-37 Mpc over this redshift range; 
the increasing dumpiness and increasing mean luminos- 
ity roughly cancel each other out. Grucially, i?max > Ri 
for all luminosities, so that isolated quasars will not be 
seriously affected by recombinations. However, the limit 
is also not extremely large, and it will clearly come into 
play for quasars that appear inside of large pre-ionized 
bubbles, preventing them from contributing as much to 
reionization as they otherwise would. Comparing to 
Figure [71 recombinations will start to limit the bubble 
growth when Xi ^ 0.6-0.8. 

5.3. The Bubble Distribution with Recombinations 

Finally, we are ready to compute the bubble size dis- 
tribution when the recombination limit is included. In 
that case, we approximate the excursion set barrier as 

This procedure makes the accumulated opacity of low-column 
density systems (with A < A^) less important, because F oc 1/R? 
along each line of sight. As a result, most of the absorption oc- 
curs near -Rmax, so the factor of two uncertainty in the mean free 
path from weak systems decreases Rmax by a significantly smaller 
amount. 
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Fig. 9. — Radius at which R = A(Ai) (i.e., the point at which 
the expected IGM mean free path equals the distance from the 
quasar), as a function of quasar luminosity. The solid and dashed 
curves are for 2 = 3 and 4, respectively. The filled triangles mark 
the luminosity-weighted mean -Rmax for each case. The dotted 
curve shows the maximum radius for an ionized bubble around an 
isolated source. 
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Fig. 10. — Bubble size distributions at 2 = 3 with and without 
recombinations (dashed and solid curves, respectively). Both sets 
assume a constant ionizing efficiency and m^nin = lOrrii. In each 
case, the curves have Xi = 0.2, 0.4, 0.6, and 0.8, from left to right. 
The filled hexagons show the fraction of space in bubbles with 
R = ^?max (see text); the Xi = 0.2 and 0.4 hexagons overlap near 
zero. 



the usual version for R < i?max- There is no barrier at 
all for R > i?max (so no such bubbles can form), and at 
R — Rmax the barrier is simply a vertical line from the 
edge of the usual version to positive infinity. Although 
this is obviously approximate, recombinations do become 
important very quickly for a given quasar luminosity, so 
it is an excellent model (but see below). With it, we can 
solve for the bubble distribution analytically (FO05). Be- 
cause there is no barrier at i? > i?max, the probability 
distribution at that size scale is a simple gaussian: 



P{S\R max 

)d6 



1 



2na„ 



exp 



2 

max 



dS, (29) 



where CTmax = f (i?„iax). The modified barrier is a vertical 
line beginning at i3(i?max)j so any trajectory that has 
(5(i?fnax) > -B(i?max) wiU be incorporated into a bubble 
at precisely the limiting radius. Their number density is 



^ rcc (^max) 



2m,, 



erfc 



S(i?max) 



V2o 



(30) 



(note that this is a true number density with units inverse 
comoving volume, not a density per unit bubble radius). 

Trajectories with (5(i?max) < ^(^max) continue their 
random walks until they cross the photon-counting bar- 
rier on smaller scales. The result must be independent of 
the trajectory at i? > i?rnax; we only care about its value 
where the barrier begins. The mass function is [excluding 
the true recombination- limited bubbles of eq. (j30p ]: 



dS p(5|-Rmax) ni,{m,z\6,R 

max J 

(31) 



where rih^m, z\6, Rmax, z) is the conditional mass func- 
tion for a trajectory that begins its random walk at the 
point (o'l'^iaxj other words, the net mass function is 



the weighted average of the conditional mass functions 
evaluated over all densities smoothed on the scale -Rmax- 

Although the integral in equation (|3ip can be solved 
analytically (FO05), the result is complicated and far 
from illuminating. We therefore simply show some ex- 
ample size distributions in Figure [TOl The dashed and 
solid curves are rih with and without recombinations, re- 
spectively. The four sets take Xi = 0.2, 0.4, 0.6, and 0.8, 
from left to right. The behavior is fairly similar to the 
analogous case of hydrogen reionization: recombinations 
are relatively unimportant for small i;, but they impose 
a sharp maximum on the bubble size distribution when 
Xi > 0.5. As expected, the limit sets in somewhat earlier 
for helium, because of its enhanced recombination rate 
and the increased dumpiness at later times. 

Unfortunately, a comparison of Figures ISland llOl shows 
that recombinations limit bubble sizes to only a few 
times the characteristic scale from stochastic reioniza- 
tion. Thus they will seriously limit the range of Xi for 
which density-driven reionization could be observable. 

There is one significant problem with our use of a single 
-Rmax value across all bubbles. During hydrogen reion- 
ization, this is an excellent approximation; each bubble 
has hundreds of sources by this stage, so fluctuations in 
the emissivity (and henc e recombination limit) are small 
(|Furlanetto et al.ll2006af ). However, we saw in ^15. 21 that 
with helium reionization we generally have only a single 
active ionizing source per bubble, and i?max is sensitive 
to its luminosity (changing by roughly a factor of two for 
each decade in luminosity). Thus some bubbles, which 
contain multiple luminous quasars or a single very lumi- 
nous quasar, will be able to grow well beyond the iimax 
shown in Figure 1101 This is simply another manifes- 
tation of our expectation that while the morphology of 
hydrogen reionization should be driven primarily by the 
clustering of ionizing sources, stochastic fluctuations are 
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much more important in helium reionization (and can 
continue to play a role late in the process). Although 
these bright sources cause only a relatively small fraction 
of ionizations, they will certainly make the recombination 
limit more of a smooth cutoff than a hard barrier. 

In principle, this formalism allows us to compute the ef- 
fective clumping factor in quasar bubbles, and hence the 
evolving recombination rate throughout reionization (see 
FO05 for a similar calculation during hydrogen reioniza- 
tion). One could then compare to the empirical ioniza- 
tion histories in ^2.1\ and constrain the importance of re- 
combinations. In practice, this is rather difficult, because 
most of the recombinations actually occur during fossil 
phases (see ^5.61 below), because we have neglected clus- 
tering of absorbers around the quasar, and because we 
lack a well- motivated model for the quasar hosts. How- 
ever, this model can be tested if one can obtain an inde- 
pendent measure of the He III ionizing background. Be- 
cause we know the emissivity from the luminosity func- 
tion, this would then constrain the mean free path - and 
hence implicitly the recombination rate and our satura- 
tion radius. 

5.4. A Consistency Check: The Post-Overlap ICM 

As mentioned above, this recombination limit controls 
the transition to the post-overlap IGM, so it is useful to 
consider observational constraints on the helium-ionizing 
photon mean free path in that regime. While our direct 
knowledge of the helium Lya forest is limited, we can 
make a number of predictions about its post-reionization 
properties from the hydrogen Lya fore st. In the optically 
thin h mit, the helium an d hydrogen iGunn fc PetersonI 
(|1965f l optical depths are (jMiralda-Escuddll993f l 

TGP.Hell _ »?thin 



The mean free path fo r a photon at the hydrogen 
Lyman-edge is therefore ()Miralda-Escudel[2003D 



TGP,HI 



(32) 



where jythin = -^HcIi/A^hi and Nj. is the column density 
of species x. Observations show a wide scatter in the 
values of r/thin in individual systems, but the me dian is 
^ 45-80 (iShull et all I^Ol iZheng etlHI l2004bD. The 
iHaardt &: Madaul (|1996f model for the ionizing back- 
ground predicts that {rjthin) ~ 40 in a fu lly-ioniz e d uni - 
verse at moderate redshifts, while Fardal et al.l ()1998D 
predict (Tythin) ~ 80. 

In the optically thick regime, a more careful calcula- 
tion including radiative transfer is required. If we de- 
fine A^LLS to be the column density of neutral hydro- 
gen for which r = 1 at the Lyman-edge of hydrogen, 
and A^HcLLS to be the column density of neutral hy- 
drogen for which r = 1 at the Lyman-edge of helium, 
the optically thin form of equation |32)) would suggest 
A'hoLLS ~ KANi^i^s/r]thin- Radiative transfer calcula- 
tions of slabs illuminated from both sides then suggest 
K ~ 2 (jHaardt fc Madaulll996f) . Thus, absorbers that are 
opaque to helium-ionizing radiation have column densi- 
ties ~ 0.1-0.2A^T,T,s. The distribution of HI absorbers 
goes like / cx N'^-^ (jPetitiean et al.l ll993). Therefore, if 
these opaque systems were entirely responsible for the ab- 
sorption, we would expect the mean free path of helium- 
ionizing photons to be Aho ~ 0.03-0. IAh. 

This is useful because the properties of 
hydrogen-ionizing p hotons are well-studied. 
IStorrie-Lombardi et al.l (1994) found that the abundance 
of Lyman-limit systems is dNi^^s/dz « 3.3[(1 -f z)/5]^/^. 



Ah « 110 



1 



Mpc 



(33) 



in comoving units, where we have again assumed / oc 
N^l^. Thus AHe 6.6, 12, and 30 Mpc at z = 4, 3, 
and 2, with about a factor of two uncertainty from the 
spread in ?7thin- These are somewhat smaller than our 
-Rmax values, but that is not surprising: even after over- 
lap, most of the universe does not lie near an active, 
bright quasar, and the mean free path will fluctuate de- 
pending on whether an active source is nearby. Other es- 
timates of the attenuation length around a ctive quasars 
after overlap are comparable to ours (e.g.. iBolton et all 
120061 take i?max = 30 Mpc). 

5.5. Helium Bubbles in the Clumpy ICM 



As described in §4.11 quasars have rather hard spectra, 
which blurs the edges of ionized bubbles (and allows some 
photons to escape to infinity). Previously we estimated 
the thickness of these bubbles assuming a uniform IGM; 
of course, a more accurate treatment includes the discrete 
absorbers discussed above. 

In detail, if the number of absorbers with column den- 
sity A^Hcii per unit redshift is /(A^hoII, -z), the effective op- 
tical depth experienced by a photon as it travels th rough 
a redshift interval (21,^2) is (|Zuo fc Phinnevlll993f ) 



Az / dAHcii/(AfHcii,^)[l-e- 



-JVHonO-HoIl(f )] 



If we then assume that /(A^hoIi) esc N ^^J^ (so that the 
helium clumps trace the H I Lya forest : IRauchll 1 998[ ) . we 
find A oc j/^'^/^, a much gentler increase with frequency 
than for a uniform IGM. In this case, half the ionizing 
photons have mean free paths smaller than « 1.9 times 
the value at the ionization edge (as opposed to a factor 
K, 3.6 for a uniform IGM). Of course, the power law 
appropriate for the Lya forest may not be applicable 
before helium reionization is complete, but it shows that 
clumping in the neutral helium phase will help sharpen 
the bubbles and hence make our model more accurate. 

5.6. Recombinations Between Quasar Episodes 

To this point, we have only discussed recombinations 
while a quasar is active. Most likely, this period fills only 
a small fraction of the age of the universe, and a substan- 
tial amount of helium can actually recombine during the 
long d orm ant phase s. We explore this aspect more fully 
in IFurlanetto et al.l (|2008( ): here we confine ourselves to 
a few words about their effect on the morphology. 

As shown in equation ([7]), the recombination time 
of mean density gas is smaller than the Hubble time 
throughout helium reionization. The recombination rate 
will be enhanced by clumping: the MHROO model pre- 
dicts that C ~ 3-4 during this era, if all the non- 
vi rialized gas is fu l ly ion ized. However, as described 
in IFurlanetto et al.l ()200^ , the clumping enhancement is 
transient during a recombination phase: once the densest 
gas has fully recombined, it can no longer aid in future 
recombinations. Thus the effective clumping for recom- 
binations quickly falls near unity. 



(34) 
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Moreover, because helium reionization occurs so 
quickly (with ~ 75% of the ionizations at z < 4), there 
is relatively little time available between quasar gener- 
ations for recombinations to occur - typically < 10% 
of th e Hubble time, except for very high redshift bub- 
bles (jFurlanetto et al.ll2008f l. Over this time span, less 
than ~ 40% of the gas is able to recombine, even with 
clumping included, and that fills only ^ 10% of the vol- 
ume. In terms of an isolated bubble, this means that 
high-density filaments and sheets will mostly recombine 
before the next quasar illuminates the region, but most 
of the volume will remain ionized, and the morphology 
will be largely unaffected - although the next generation 
of quasars will have to reionize these dense pockets. 

On the other hand, ionized regions generated by rare, 
extremely high-redshift quasars will recombine signifi- 
cantly before being illuminated again; as a consequence, 
many of their photons are essentially wasted during this 
early phase. 

6. DISCUSSION 

We have examined several aspects of helium reioniza- 
tion using_sini£kanal^^ models. We first showed that 
the lHopkins et al.) ()2007d( ) luminosity function, together 
with standard template quasar spectra, produces ~ 2.5 
ionizing photons per helium atom at z = 3, nicely consis- 
tent with observational hints that the IGM truly becomes 
fully ionized at that time. Indeed, even with an average 
clumping factor C ~ 3, the known quasar population 
can ionize helium by that time. Interestingly, the quasar 
emissivity increases rapidly at z > 3, with ~ 75% of the 
ionizing photons produced after z = 4. This is somewhat 
faster than naive models based on the collapse fraction 
predict, implying that supermassive black holes in quite 
massive halos dominate the photon budget. 

In our calculations, we have used a quasar luminos- 
ity function with a shallow faint - end s l ope at z > 3 
(iHunl. et al.l l2004l : ICristiani et"all 120041 : iHopkins et all 
l2007df) . This implies that most ionizing photons are pro- 
duced by bright quasars {L > L*), which are of course 
relatively rare. As a result, the initial phases of reion- 
ization are dominated by rare, large bubbles (with char- 
acteristic size i?char ~ 15 Mpc) whose overlap is ran- 
dom. We described this regime with a "stochastic" model 
where sources are distributed according to Poisson statis- 
tics. In this model, the bubble size distribution follows 
a power law (roughly), so there is no true characteris- 
tic scale - although the power law is so steep that, in 
practice, nearly all bubbles contain one or at most a few 
sources anyway. Thus, in the early stages of reioniza- 
tion, we expect bubbles to have a characteristic size cor- 
responding to a single or at most a few sources. How- 
ever, more recent data shows a somewhat steeper faint- 
end slope (iFontanot et al.l 120071 : iBongiorno et al.ll2007t 
ISiana et al.ll2007t ). This leads to a broader distribution 
of luminosities, making the stochastic phase somewhat 
less important. 

Large-scale clustering of the ionizing sources will even- 
tually cause the bubbles to grow beyond this point. We 
use the model of FZH04 to describe this phase, in which 
large-scale overdensities are ionized first. This imprints 
a well-defined characteristic scale, which first exceeds 
the size of the stochastic bubbles when Xi 0.5 and 
would exceed 100 Mpc by the end of reionization if al- 



lowed to grow unchecked. The scale is nearly indepen- 
dent of redshift at a fixed Xi, but it does depend on the 
clustering of the quasar hosts, with more massive host 
galaxies driving larger bubbles. In our calculations, we 
used a simple halo-based model to represent quasar hosts. 
This is reasonable if the observed relation s between black 
hole mass and halo properties ( Magorr ian et al.l 119981 : 
iGebhardt eralll2000l : iFerrarese fc Merrittll2000(n iold at 
high redshift as well. But many detailed models of the 
quasar luminosity function rely on mergers to trigger 
quasar phases. This would change the distribution of 
quasar luminosities over halos an d so change some of th e 
predictions of our model (see also lCohn fc Chaiigll2007| ). 
We intend t o examine these effects w i th the semi-numeric 
approach of iMesinger fc Furlanettol ()2007[ ). 

However, even during this phase nearly all quasars 
exist in isolation: provided that their duty cycle 
tQSoH{z) ^ 1, very few bubbles will contain more than 
one active source. Thus the ionizing background around 
each source is easy to compute (modulo the nearly uni- 
form background of high-energy photons with large mean 
free paths). Together with a model for the IGM, this al- 
lows us to estimate the mean free path of helium-ionizing 
photons as they propagate through the ionized bubbles 
and specifically to compute how far they can travel be- 
fore being consumed by dense, optically thick absorbers. 
We find that, at the ionization edge, -Rmax ~ 35-40 Mpc 
when averaged over the entire quasar population. Once 
bubbles grow to this size, only a fraction of the photons 
would ionize new helium atoms. As described in FO05, 
the effective clumping factor of the IGM will increase 
as more and more photons get consumed by these dense 
systems, so reionization will slow down. Eventually, the 
ionized gas distribution will map onto the more uni- 
form background characteristic of the post-reionization 
era (although ev en then large fluctuations will persist; 
iShull et al.li2004J) . 

Our model does not include recombinations within 
bubbles that contain zero active quasars over a period 
of time, except as an overall_ scaling factor to the the 
ionizing efficiency (through A^rec)- The helium recom- 
bination time is short enough, and the IGM dumpi- 
ness large enough, that a significant fraction of helium 
atoms (mostly in high density regions) can recombine, 
although > 90% of the volume will remain highly ion- 
ized. We examine thes e "fossil" bubbles more closely in 
iFurlanetto et all ()2008f ). 

Compared to hydrogen reionization, the helium era 
passes through the same qualitative phases: stochastic, 
isolated bubbles followed by a rapid increase in the bub- 
ble size when clustering becomes significant and finally 
a "saturation" when the mean free path of an ionizing 
photon becomes smaller than the bubble size. However, 
quasars are much more rare than star-forming galaxies 
at z > 6 (at least in most models of the high-redshift uni- 
verse) , so the early stochastic phase is much more impor- 
tant; for hydrogen, it ends by Xi ^ 0.1. There is thus a 
much shorter window in which density-driven clustering 
dominates the morphology, when -Rstoch < R < Rmax- 

Because there is typically no more than one quasar 
per bubble, the scale i?max is much easier to estimate 
than for hydrogen. FO05 estimated it during hydrogen 
reionization by balancing the total recombination rate 
within the bubble against the emissivity of the ionizing 
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sources. This depends on the recombination history and 
uncertain IGM density distribution, and it assumes that, 
within each bubble, the ionization always begins at low 
densities. As a result, i?max is quite sensitive to many of 
the input parameters (varying by about a factor of two 
between case- A and case-B recombination, for example). 

For helium reionization, i?max can be much more ro- 
bustly determined because it depends only on ionization 
balance around a source population (whose luminosity 
function is quite well-known, at least at z < 4) in an IGM 
whose density structure is extremely well-determined by 
the H I Lya forest (although the spread in luminosities 
also makes it more difficult to incorporate directly into 
the simple analytic model). The helium reionization era 
will thus offer a much more stringent test of the tran- 
sition to the recombination-dominated phase than the 
hydrogen era. 

There is, however, a possibility that hydrogen and he- 
lium reionization are more similar than it may appear. 
We showed that observations at least allow the possibil- 
ity that 2: w 3 galaxies produce nearly as many helium- 
ionizing photons as quasars. Their importance for reion- 
ization depends entirely on their escape fraction, which 
is difficult to estimate. If they do contribute, we would 
expect a slower evolution of Xi with redshift, a softer 
ionizing background, many more small ionized bubbles, 
and a more uniform, smaller mean free path (and hence 
-Rmax) for the helium-ionizing photons. 

We have focused exclusively on analytic models here, 
but of course a full picture of helium reionization requires 
numerical simulations that capture both the stochas- 
tic and large-scale clustering elements of our model, as 
well as the complex geometry of the IGM and source 
distribution. The models described here can be ex- 
tended to "hybrid" semi-numeric schemes ()Zahn et al.l 
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